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Q_ij The rate of monopole decay into a dyon and an electron in an inhomoge- 

O 1 neous external electric field is calculated by semiclassical methods. Comparison 

is made to an earlier result where this quantity was calculated for a constant 
, T"! \ field. Experimental and cosmological tests are suggested. 
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1 Spontaneous and Induced Decay in Exter- 
nal Fields 



Spontaneous non-perturbative processes of particle production in QFT 
(Schwinger processes), or vacuum decay processes have been studied since the 
historic papers of Euler-Heisenberg and Schwinger [1, 2] on e + e~ generation in 
a constant electromagnetic field. Voloshin, Kobzarev, Okun [3] were the first to 
treat false vacuum decay in a scalar field theory with a stable and a metastable 
vacuum states. Later Callan and Coleman [4, 5] gave this problem a 1-loop 
treatment, calculating both the exponent and the preexponential factor of false 
vacuum decay probability. 

Consideration of induced [6, 7] non-perturbative particle creation was a natu- 
ral extension of the scope of the problems described above. The term "induced" 
denotes the situation in which the initial state is not vacuum, but rather contains 
some particle(s). False vacuum decay in a scalar and spinor field theory, induced 
by presence of an external particle acting as a "catalyst" or "nucleation center" , 
was treated semiclassically up to 1-loop preexponential in [8]. 

On the other hand, generalization of Schwinger processes description can be 
thought of as extending the class of fields in which the appropriate process takes 
place. The original Euler and Heisenberg calculation in QED was performed for 
a constant field. For harmonic plane waves calculations had first been done by 
Schwinger in the cited paper. One can make sure [9] that the same expression 
is true for adiabatically varying fields. Narozhny and Nikishov [10] calculat- 
ed exactly the effective Lagrangian in an electric field, dependent on time as 
E(t) ~ cos ^(m) • ^ semiclassical treatment of a broad class of fields was given 
in [11]. Semiclassical methods were further developed basing upon WKB ap- 
proximation [12] and the so called "worldline instanton method" [13, 14, 15]. To 
mention some other exact results, Fried and Woodard [16] gave an expression 
for arbitrary light-cone coordinate dependent field E(x$ ± xi). For a compre- 
hensive review of recent developments in Euler — Heisenberg effective actions 
the reader may consult Dunne's review [17]. This extensive list (most part of 
which has been left behind in order not to overload the reader) of Schwinger 
processes in inhomogeneous external fields is mostly related to QED processes. 
Some authors also dealt Hawking radiation in a Schwinger-like manner [18, 19]. 
It should be emphasized that most of the papers on inhomogeneous field vacuum 
decay consider spontaneous processes. 

In the present short paper we suggest combining the both generalizations of 
Schwinger processes. The possibility of an induced monopole decay into a dyon 
and an electron was first suggested in [20]. Later it was calculated in a constant 
electric field in [21] up to leading classical exponential factor. This problem 
deserves attention per se, but below we also give reasons for astrophysicists to 
be interested in such processes. 

The paper is organized as follows. In section 2 we describe the general tech- 
niques of dealing particle decay in terms of Euclidean 1-particle path integral. 
The sub-barrier trajectories and the leading exponential factor are calculated in 
section 3. The circumstances under which the process considered might become 
significant for observers, are investigated in section 4. In section 5 we discuss our 
results. 
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Figure 1: Full 1-loop Green function of a monopole is obtained by summing over 
the insertions of electron-dyon loop into monopole's Euclidean trajectory. 

2 Quasi-Classical Approximation to Path 
Integral 

We are going to study 't Hooft-Polyakov (non point-like) monopole and dyon. 
The masses of these particles M m and M<j are of order of the scale where M\y 
is generally the scale at which spontaneous gauge symmetry breaking takes place 
(the mass of W-boson), a coupling constant. At the same time, their sizes are of 
the order of magnitude , thus in weak coupling limit de Broglie's wavelengths 
of monopole and dyon are far smaller then corresponding sizes, therefore, these 
particles are essentially classical objects. 

On the other hand, 't Hooft-Polyakov monopole does not possess a well- 
established local field-theoretical description. That is why it is reasonable to 
treat monopole and dyon in terms of 1-particle theory (quantum mechanics), 
evaluating (bosonic) Feynman path integrals semiclassically with restriction on 
the trajectories of classical motion r 3> , where r is a typical trajectory size. 
We are going to study monopole propagator in imaginary (Euclidean) time. To 
do that semiclassically, one should first find closed-loop trajectories in Euclidean 
time, and then calculate determinants, corresponding to quantum fluctuations 
around them. It is essential to take into account only closed loops contributions 
because only paths with finite classical action on them are relevant. 

To find the monopole decay probability in an inhomogeneous external field 
one has to calculate corrections to its propagator in the presence of external 
electric field. What kind of trajectories will those corrections correspond to at 
the classical level? As monopole is now not forbidden (due to the presence of the 
field) to decay into an electron and a dyon, these are configurations containing, 
beside the monopole pieces, dyon and electron pieces. Therefore, finding the 
full Green function of a monopole in an external field is equivalent to evaluating 
Green function with arbitrary number of electron-dyon loop insertions 

G(T, 0) = (T, 0) + G* 1 ) (T, 0) + (T, 0) H 

The free propagator G^ (T, 0) corresponds to the diagram without dyon-electron 
loop insertion, G^ (T, 0) denotes Green function with one insertion of electron- 
dyon loop etc., see Fig. (1). Note that the diagrams of Fig.l are not Feynman 
graphs, but are simply the classical trajectories in OX4X3 plane. However, re- 
summation of all electron-dyon loop insertions resembles closely the analogous 



resummation of all one-particle irreducible vacuum polarization diagram inser- 
tions. 

Taking into account that the electric charge of the magnetic monopole is zero 
and using the first quantized approach, the propagator of the monopole in the 
Euclidean time can be written as follows 




where y(t) is the monopole variable. The first correction comes from Euclidean- 
time configuration, shown in Fig. 2. Note that electron and dyon can go round 
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Figure 2: One-loop Euclidean configuration ("instanton" with winding number n = 1). 

the loop multiply, winding over it with some respective winding numbers m, n. 
This corresponds to the dashed trajectories in Fig. 2. Closed trajectories 
classified according to their winding number, will be called, following Dunne 
and Schubert [14], "worldline instantons". They are analogous to Yang-Mills 
instanton-antiinstanton pairs. Propagator correction becomes after resumma- 
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Figure 3: One-loop Euclidean configuration with multiply winding trajectories ( "world- 
line instantons"). Dashed bold arcs correspond to extra winding paths compared to 
Fig.2, arrows indicate winding direction. 

tion over winding number m, n 

m,n 

Here x = (0,0,0,0), z = (T, 0,0, 0) are monopole initial and final Euclidean 
coordinates. Equivalent notation G(x, z) = G(T, 0) will be used for such config- 
urations. 



To give the reader the feeling of how S^ 1 '^ can be organized, spontaneous 
vacuum decay case in QED may be brought as an example. Configurations with 

(n) 

greater winding number n would give, in general, S d ~ nS c i, i.e. n times greater 
classical action, and the preexponential factor would depend on n as K n = ^jK\. 

In the present case, it is reasonable to expect S^ 1 '^ ~ mS^ + nS^ . 

This resummation is not being performed here, only first term K\^e ci 
in the sum being left. Below K± t i will be denoted simply K, S± t i as Sd- This 
resummation would be necessary only in case of S c i — > 0, otherwise the contri- 
bution of higher winding number terms is strongly suppressed. As S c i — ► 0, both 
the approximation of point-like monopole and the semiclassical approximation 
break down, so the interest to this resummation is mostly hypothetical. 

For the lowest winding numbers m = 1, n = 1, the first correction due to 
electron-dyon loop is given as 



!) (T, 0) = J Vye~ Mm S ^ dT VxVze~ s ^ z ^ , (2) 

where action S[x(u), z(v),A] is the common action for the charged particles (dyon 
and electron) in the external field A^, x(t) and z(t) are respectively electron and 
dyon coordinates, their worldlines being parametrized by u, v 

S[x,z,A] = m j'Vx^du + ie J A ext (x)xdu + M d j Vl^dv - ie J A ext (z)zdv . 

Monopole coordinate is easily integrated out, as monopole does not interact with 
the electric field, 

j^^^-^j^.^^^, (3 ) 

here h is the vertical size of the loop (see Fig. 2). 

The remaining path integral has a more complicated structure of modes. It 
always has zero modes, corresponding to the shifts of the loop. This will con- 
tribute a preexponential factor of J[x(t), z(t), A], which is the Jacobian, arising 
due to the substitution of variables: normalized zero modes coefficients to the 
collective coordinates (position of the loop). 

The other part of the preexponential factor will come from the non-zero 
modes. Among them there will be at least one negative mode, corresponding 
to the overall inflation of electron-dyon loop. Presence of extra negative modes 
is not evident without a special investigation [21]. In general, the path integral 
yields 

G« (T, 0) = J d 4 yG {0) (x, y)G<® (y + Ay, z)Ke' s ^ (4) 

where S c i is the action of dyon and electron in the external field on a classical path 
with proper boundary conditions (see Fig. 2). The integral over d 4 y emerges due 
to integrating over all positions of the loop, and the two free Green's functions 
belong to the purely monopole trajectories. Shift Ay accounts for non-zero size 
of the loop, and can be neglected as the condition T 3> Ay is imposed. 
K contains contributions from the Jacobian and from non-zero modes 



K = J[x(t),z(t),A] 



detF \ 
detF ) 



Fluctuation operator F can be found in [21], whereas Fq is normalization oper- 
ator, analogous to operator of fluctuations around trivial vacuum configuration 
in false vacuum decay problems. 

We have just given a sketch of structure of the whole expression for the prop- 
agator, as below we are going to be interested mainly in the leading semiclassical 
exponent, leaving the preexponential part for a more sophisticated analysis in 
future. 

How can (4) be put into direct correspondence to the mass shift of the 
particle? If the full Green's function is G{x,z) « t^jj / tf+rrfi+Sm'* ' mass sn ^ 
Sm 2 being imaginary or real, then expanding in the powers of 5m 2 one gets for 
the variation of Green's function 

SG(x, z) = G {1) (x, z) = -5m 2 j d 4 yG (0) (x, y)G (0) (y + Ay, z) 
Comparing this with (4), one makes sure that 

5m 2 = Ke~ Scl 



3 Exponential factor 

We are going to do path integral over electron and dyon coordinates by steepest 
descent approximation, applying the ideas of world-line instanton method due to 
Dunne and Schubert [14, 15] to induced decay problem. To find the exponential 
factor of the probability one should solve classical equations of motion for dyon 
and electron in an external electromagnetic field, find closed configuration of 
Euclidean trajectories and minimize the action on them with regard to trajectory 
parameters. A single-pulse electric field directed along OX3 axis 

^3 = Eq — 2 

cosh (cot) 

will be considered. 

Equations of motion will be of the form 

m T~~7% = -ieFu, v (x)x„ 

M d ^^==ieF^ u {z)z v 

where m, M d are electron and dyon masses respectively. Here for convenience, 
dimensionless parameters 7 = a = jj-, (3 = are introduced. Solution to 

these equations for the field ^3(2:4) = — tan(c<JX4) will be for the electron and 
for the dyon respectively 

m 1 

X3{u) = — — arcsinh (7 cos(2-7ru)) — a 

e ^7\A + 7 2 

m 1 . / 7 \ 

£4 (it) = — — — arcsm — sm(27ra) , it 6 [— ito,ito] 

eE i \V 1 + 7 2 J 

Tfl 1 / 7 \ 

z^iv) = , arcsinh I — cos(27ri;) I + b 

eE L . 7 2 \a / 

Tfi 1 / 'y \ 

za{v) = — -arcsin - sin(27ri>) ) ,v € [-v ,v ] , 

e^7 \y/a 2 + 7 2 / 
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Figure 4: Trajectories in the Euclidean plane (£3, £4) for various values of parameters 
0, 7; a = 5 • 1CF 3 , taken as length unit. 



where a, 6, iio, t>o are some constants. If we calculate the action for these trajec- 
tories we will receive a four-parametric expression, depending on a, b, uq and vq. 
We should find the values of these parameters which minimize the action. One 
can minimize the action using the standard procedure and taking into account 
that resulting loop should be closed, i.e. x^(uo) = ^4(^0). Boundary conditions 
will in this case enter the equations as some constraints with Lagrange multipli- 
ers. Instead of this we can use mechanical analogy (see [7]), for which m, M d and 
M m are forces acting in the vertex, so, the minimum of the action corresponds to 
the value of parameters for which we have equilibrium in the vertex. Condition 
of the equilibrium is 

m cos a e (uo) — cos Od(^o) = 0, 

msma e (u ) + M d sma d (v ) = M m . (5) 

a e (uo), ad(vo) are the inclination angles of corresponding trajectory in the vertex. 
For the electron we have 

tano e = \/l + 7 2 cot 2ttuq. (6) 

Using analogous relation for dyon and (5) one can easily find 

, 2 /, , 2 , ((M m + mf - Mf) (Mi - (M m - mf) 
tan z 27rii = (1 + 7 )- — -9 -, 

f\ {(M m + mf - Ml) [M\ - (M m - mf) 
«V (Ml + Ml-m 2 ) 2 

It follows immediately from here that for the boundary conditions to be properly 
satisfied, dyon mass should lie within the interval M m — m < M d < M m + m. 
The action on the extremal trajectory is 



tan 2ttvq 



Sd(l) = S e (j) + ^(7) - 2M m Xi(u ), 



(7) 



where 
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1 27rm 2 
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[2n + 



[2^0 + g] 



where square brackets [x] in the last expression denote minimal integer part of x. 
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For sina e (i > 0, or tan27ru > and tan27rv > (see (6)) i.e. Mj < m 2 + 



it follows that 
1 



u 



2vr 
1 

2^ 



arctan 



(1 + 7 2 
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1 + 
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so, the action reads as 

1 2vrm 2 

5 d (7) = 



7 2 eE 

1 2ttto 2 
T 2 e £ 



2«o \/l + 7 2 ~~ — arctan ( tan j 



2/ o\/ 1 + - — arclau 



1 

7T 



tan 2ttvo 



(8) 



For the case Mj>m 2 + it follows that 



1 1 

un = arctan 

2 2vr 



(I + 7 2 



((M m + m) 2 - Mj) (Mj - (M m - m) 2 ) 1 



(mi + Ml-Mj)' 



Corresponding value of the action is 
1 2vrm 2 



Se(l) 



7 2 eE 



1 + 7 2 arctan 



7T 



tan 2-kuq \ 



(9) 



A family of extremal trajectories for several values of 7 is shown in Fig. 4. One 
can see the overall plot of S(^) versus Keldysh parameter 7 in Fig. 5. 

The relation jH, where To is monopole mass imaginary part (width) in a 
constant electric field (7 = 0) can be plotted and is shown in Fig.6. One can 
easily see the increase of the particle width with the increase of 7. This increase 
becomes stronger when the difference between M^ and M^ grows. 

Action S tends to zero as 7 — ► 00, therefore, at some point, resummation of 
higher winding number contributions (1) will be necessary. At sufficiently low 
S the whole process should rather be dealt within common perturbation theory, 
because semiclassical approximation to the non-perturbative mass shift becomes 
invalid here. 

As it was mentioned in section 2 one has a criterion for applicability of the 
first quantized approach to dyon and monopole. Namely, the loop should be 




Figure 6: ^- dependence on Keldysh parameter 7 for temporally inhomogeneous field. 



large enough for dyon to be treated as a point-like object. That is, characteristic 
loop size, which is about the electronic loop "radius" 

-v fTE' 7>>1 



7 eE' 



7 << 1 



should satisfy the condition 
that is, 



r e > M w \ 



oo < M w 

E ^ ^ M w 



Here £" cr = ^- is the critical field value for electron. For imaginable processes 
in cosmology, field switch-on rate is obviously less then the enormous 100 GeV 
~ 10 26 Hz, thus in fact the first criterion is always satisfied. The second criterion 
means that E « 10 8 .E cr , which is valid for most of magnetic stars and Reissner- 
Nordstrom black holes (see section 4). 

However, from general arguments it becomes clear [14] that semiclassical 
approximation is wrong when S <C 1. Therefore, this general limit could give 
us a more strict criterion of applicability of the formula (7) than monopole size 
considerations. 



3.1 Preexponential: Negative modes 

Here we limit ourselves to qualitative considerations. The presence of the zero 
mode was already discussed and used in the resummation above. We might 
state, due to the fact that for 7^0 one has a situation, totally analogous to 
that described in [21]. At sufficiently small 7 if Mj < M 2 L + m 2 one always has 
one negative dilatational mode, corresponding to overall inflation of the loop. It 
provides a possibility for monopole mass to acquire an imaginary part. If M 2 > 
+ m 2 , another negative mode comes into existence (see [21]), thus making 
loop contribution to the monopole mass real and describing mass renormalization 
of a stable monopole. In future, we are going to elaborate a numerical criterion 
for 7, below which the above statements are true. 



3.2 Spatially inhomogeneous field 

Dunne and Schubert have shown [14] that in a spatially inhomogeneous field 
particle decay probability can be obtained by an analytic continuation of the 
result for the field with the same temporal inhomogeneity 

7 — ► 27. (10) 

See also [31, 32] for a more detailed discussion of this analytic continuation 
procedure. The general fact on spontaneous (induced) processes in spatially in- 
homogeneous fields is that the imaginary part of the effective action (particle 
width) decreases with the increase of 7, becoming zero at some point. Physically 
one can interpret this fact easily: the field characteristic size becoming too nar- 
row, so that no pair of particles can travel a path long enough to gain energy 
necessary for leaving the barrier. The explicit formula for the action reads 



Sdil) = 5 e (»7) + $$7) - 2M m i 4 («o,«7). 



(11) 



where 

r, r ■ \ 1 27rm 2 / /- =■ 1 / tan 27ruo \ \ 

SeXfn) = -2 — W -2-uoVl-T ^ + -arctan 

T ^ \ 7T \ Vi-7 2 / / 

„ . . , 1 2vrm 2 / „ / j* 1 / tan 2irv \ \ 

Following the simple analytic continuation rule, we show how the exponential 
factor behaves for field 
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Figure 7: Action S of electron-dyon loop in units of ^ versus Keldysh parameter 7 
for spatially inhomogeneous field. 

In Fig. 7 one can see that the process of monopole decay becomes infinitely 
suppressed as 7 — > 1. At higher values of inhomogeneity parameter 7 > 1 decay 
is forbidden. This corresponds fully to the earlier results of [14] for spontaneous 
pair creation. 

Below in Fig. 8 we plot the ratio of monopole width in a spatially inhomoge- 
neous field with parameter 7 to its width in the constant field case. 




Figure 8: j£- dependence on Keldysh parameter 7 for spatially inhomogeneous field. 



4 Monopoles in Cosmology and High- 
Energy Physics 



The issue of magnetic monopoles, present in almost every 4-dimensional non- 
abelian gauge theory with a scalar Higgs field [22, 23] has long been an important 
problem in high-energy physics. Violation of baryon number conservation law 
in presence of a magnetic monopole [24], charge quantization [25] are just a few 
well-known features of rich monopole physics. 

One of the problems of monopole physics is the overestimation of monopole 
concentration in the Universe with regard to experimental limit. E.g. it has 
been pointed out that our Universe should contain relatively large concentration 
of free primordial Polyakov - 't Hooft monopoles of M ~ My/ /a ~ 10 4 GeV (over 
10 orders of magnitude larger than the upper experimental bound), unless they 
are bound in meson-like states or there are alternative mechanisms of monopole 
decay. GUT monopoles' (M ~ 10 17 GeV) generation just after the spontaneous 
violation of the gauge symmetry of the GUT group was long being thought to be 
a major trouble for the standard cosmological model, until having been solved 
in terms of inflationary cosmology in [26] . 

"Phenomenology" of monopoles nowadays has set some limits on their abun- 
dance. The observed flux J on Earth is limited from above as J < 1.4 x 
1(T 16 ^ — [27]. The so-called Parker bound coming from magnetohydrody- 
namic considerations limits the flux of monopoles by J < 1 x 10" 15 — ^ . 

1 J cin z sec sr 

The expected velocities of the monopoles lie in the range 10 -4 < v/c < 1, masses 
of "primordial" GUT monopoles are expected to be of the order of magnitude 
10 17 GeV or higher, monopoles with masses 10 7 . . . 10 15 GeV could have been 
generated during the later phase transitions (after breaking the GUT symme- 
try). Lowest 't Hooft monopole mass is reached in BPS limit M ~ My/ /a. The 
limit on concentration of monopoles depends on their velocity and flux, and can 
roughly be estimated tim < 10~ 26 • • • 10~ 22 cvrT 2 ' . 

Thus the concentration of monopoles in the Universe is a very important 
quantity, as we have experimental bounds on it and, on the other hand, the 
presence of monopoles may have influenced its evolution. Therefore, it is very 
important to know well the possible channels of monopole decay. In particular, 
one is interested in external field induced decay processes due to several reasons. 

First, one cannot fully eliminate the scenario of monopoles decaying within 
stellar matter suggested long ago by Zeldovich and Khlopov [28]. One should 
pay particular attention to these processes in extremal Reissner-Nordstrom black 
holes, as these objects may create electric fields close to the critical value. (The 
upper boundary of Reissner-Nordstrom black hole electric field is actually limited 
from above by 2 x 10 12 Gs [29]). Here one can use constant field limit 7 = 0, 
and the value of the typical action (7) is S ~ 10, so that monopole decay is not 
infinitely suppressed by the factor of e~ s ~ 10~ 4 ' 3 . 

Second, at the start of inflation stage when monopoles had just been born, 
we may have some non-zero expectation value of (electromagnetic) field, which 
might have catalyzed monopole decay into something else. 

Third, one could consider decay of monopoles in the interstellar media, as it 
has been argued that magnetic fields of the order of magnitude 10~ 2 Gs must 
exist in dense molecular regions and IR nebulae [30]. 

Let us estimate realistic values of Keldysh parameter 



1(T 19 sec E cr 

for possible applications of our technique, where r is characteristic inhomogeneity 
time, E typical field value, and E cr critical field value for electron (4 • 10 13 Gs). 
The most rapid processes observed in the modern Universe have to do with 
pulsars and can have r ~ 10 -3 sec. On the other hand, typical magnetic field of 
a pulsar can reach the order of 10 _1 S cr . Thus 7 is extremely small and one can 
use the stationary approximation from [21]. 

In terrestrial conditions, lasers with r ~ 10~ 15 • • • 10~ 16 , and E ~ 
10~ 3 • • • 10~ 5 could be within the reach of modern experimentalists. Typical 
action in this case is S ~ 10 3 , which still suppresses monopole decay. If gamma- 
lasers could produce pulses short enough, one could hope for diminishing r by 
2 — 3 orders of magnitude and reach values S ~ 1. However, such parameter- 
s are out of reach at present time. Moreover, the semiclassical approximation 
breaks down at such values of S, as we have mentioned already. Speculating 
further, one can put a little fantasy in it and conjecture that monopoles could 
be produced at some high energy facility, and then directed into a device, where 
the corresponding time-dependent field is created. Then monopole decay into a 
dyon and an electron could be observed. However, we understand both theoret- 
ical and experimental difficulties of implementing such a project. All attempts 
to detect monopoles at accelerators have failed up to now. E.g., production of 
(Dirac) monopoles at Tevatron has been considered in [33]. From the data on 
the existing facilities, bounds upon monopole mass and cross-section have been 
established: M > 300 GeV, a < 10~ 37 cm 2 . This imposes bounds upon the pos- 
sible monopole production at, e.g., LHC. If its luminosity be ~ 10 34 cm _2 s _1 , 
the upper limit of monopole production would be set as 10 4 per year, if we rely 
on the upper bound for cross-section. It is obviously a non-trivial task to detect 
monopole production itself, and even more complicate done to realize its decay 
under the conditions assumed in this paper. 

5 Discussion 

In this short Note the description of the induced monopole decay has been gen- 
eralized to the non-stationary field case. Comparing to stationary field config- 
uration, one can conclude that, as expected, monopole decay is enhanced by 
a temporally-inhomogeneous field and suppressed by a spatially inhomogeneous 
field. It has been shown that, despite being non-perturbatively suppressed, this 
process may take place under some exotic conditions, e.g. in Reissner-Nordstrom 
black holes and in pulse gamma or x-ray lasers. 

Non-stationarity of field becomes a key factor in the latter case, allowing the 
classical action on the Euclidean closed worldline to become sufficiently low and 
thus cease to suppress monopole decay. On the other hand, we have shown that 
even for the most rapid processes in cosmology, it is generally possible to use 
constant field approximation from [21]. 
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